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Abstract 

The thermodynamic stability of onion-like carbon (OLC) nanostructures with respect to highly oriented pyrolytic graphite 
(HOPG) was determined in the interval 765-1030 K by the electromotive force {emf) measurements of solid electrolyte 
galvanic cell: (Low) Pt|Cr3C2,CrF2,OLC|CaF2S.c.|Cr 3 C2,CrF2,HOPG|Pt (High). The free energy change of transformation 
HOPG = OLC was found positive below 920.6 K crossing the zero value at this temperature. Its trend with temperature was 

well described by a 3 rd degree polynomial. The unexpected too high values of ( — ' T ) =Ac P (T) jointly to the HR-TEM, 



p 



STEM and EELS evidences that showed OLC completely embedded in rigid cages made of a Cr 3 C 2 /CrF 2 matrix, suggested 
that carbon in the electrodes experienced different internal pressures. This was confirmed by the evaluation under constant 

volume of — by the - ratio for OLC (0.5 MPa K -1 ) and HOPG (8 Pa K -1 ) where a and k are the isobaric thermal expansion 

dT k 

and isothermal compressibility coefficients, respectively. The temperature dependency of the pressure was derived and 
utilized to calculate the enthalpy and entropy changes as function of temperature and pressure. The highest value of the 
internal pressure experienced by OLC was calculated to be about 7 GPa at the highest temperature. At 920.6 K, A,H and 
A r S values are 95.8 kJ mol -1 and 104.1 JK _1 mol -1 , respectively. The surface contributions to the energetic of the system 
were evaluated and they were found negligible compared with the bulk terms. As a consequence of the high internal 
pressure, the values of the enthalpy and entropy changes were mainly attributed to the formation of carbon defects in OLC 
considered as multishell fullerenes. The change of the carbon defect fraction is reported as a function of temperature. 
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Introduction 

The modifications of the internal arrangements and related 
energies of single wall carbon nanotubes (SWCNTs) in bundles 
were studied [1] by our group in the particular situation where 
their dilatation due to the high temperatures was hindered being 
the bundles embedded in a matrix much less dilatable. This 
generates high internal pressures, which are a direct consequence 
of the increase of the internal energy of the system producing 
changes in the SWCNT configuration inside the bundle and 
deformation of SWCNTs too. 

The knowledge of the thermodynamic stability of nanostruc- 
tures is a fundamental aspect to design new nanostructured 
materials as well as to forecast their behaviour. Most of the 
information on this subject comes from computational studies 
[2] [3] [4] [5] [6] and to a lesser extent by experimental works 
[7] [8] [9] . To the best of author knowledge and with the exception 
of reference [10] [11] [12] there are no articles dealing with 



thermodynamic measurements at high temperature on carbon 
nanostructures. 

On the other hand, the 3D carbon phase diagram, where the 
third axis is the particle size, is still under consideration. The 
contribution of the surface energy to the bulk Gibbs free energy 
per atom of a cluster of n atoms in a given phase is strongly 
affected by the number of atoms itself. In general, a review of the 
literature regarding the structure of carbon nanoparticles [13] 
highlights that at sizes below 1.8 nm, other carbon forms are 
abundant, such as fullerenes and onion-like carbon (OLC). 
Therefore, Kuznetsov et al. [14] suggested to assign a corre- 
sponding region of the phase diagram to closed-shell ,sjf) 2 -bonded 
nanocarbons. In their phase diagram the regions of stability are 
specifically indicated for fullerenes and OLC for n = 1 0 to 1 0 2 
atoms. The complete and thorough vision of the phase diagrams of 
these nanostructures in comparison with the phase diagram of the 
corresponding bulk materials has not been realized yet. There is 
still a great deal of work to be done in the construction of a 
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complete P-T—n phase diagram of carbon. It is likely however that 
knowledge gained from theoretical and computational studies [15] 
of nanocarbon stability plays an important role in the conception 
of such a diagram, with thermodynamic treatments of phase 
equilibrium showing the right way for pursuing the final goal. 

Onion-like carbon are another class of carbon nanostructures, 
which are not exhaustively investigated yet though some published 
results envisage interesting applications such as an easy route to 
produce diamond nanoparticles [7] [8] [16] as well as to use them 
as reaction cells on the nanoscale [9]. 

The scope of this work is to study, in conditions of reversibility, 
the thermodynamics of the high temperature transformation 
under volume constraint from highly oriented pyrolytic graphite 
(HOPG) to OLC. 

Experimental 

Since the experimental set-up was published elsewhere 
[10] [11] [12] and there is the will of the authors to give in this 
paper more room to the chemical thermodynamics of our 
experiment and discussion of the results, we decided to describe 
more in details the wide experimental part in the Supporting 
Information S 1 . A list of all the aspects that took part to the whole 
experimentation of this work, which are reported in the 
Supporting Information SI, follows: 

1. Starting materials for the electrode preparation 

1.1. HOPG, OLC, Cr 3 C 2 , CrF 2 

1.2. Characterization 

1.2.1. Thermogravimetry (TG) and Differential Thermal 
Analysis (DTA) 

1.2.2. x-ray diffraction (XRD) 

1.2.3. micro Raman (mR) 

1.2.4. High Resolution Transmission Electron Spectroscopy 
(HR-TEM) 

1.2.5. Electron Energy Loss Spectroscopy (EELS) 

1.2.6. X-ray Photoelectron Spectroscopy (XPS) 

1.3. Preparation of electrodes and cell assembly 

1.3.1. Electrodes 

1.3.2. Cell assembly 

1.4. Experimental apparatus for the emf measurements 

1.5. Procedure adopted for the emf measurements and their 
data acquisition 

1.6. Analysis of electrodes before and after experiment 

1.6.1. X-ray diffraction (XRD) 

1.6.2. microRaman spectroscopy (mR) 

1.6.3. EELS 

Results and Discussion 

In order to satisfy the requirement of reversibility, the 
experimental technique utilized is the electromotive force (emf) 
measurement as function of temperature of galvanic cells, with 
CaF 2 single crystal as solid electrolyte. Since long time this 
technique was utilized in our laboratory and it is well known that 
the method is one of the best way to achieve reliable 
thermodynamic data. 

The galvanic cell: 

(Low) Pt | Cr. j C 2 ,CrF 2 ,OLC | CaF 2 s.c. | Cr,C 2 ,CrF 2 ,HOPG | Pt 
(High) (A) 

was assembled to perform the emf vs T measurement. Cell A is a 
fluorine concentration cell being the solid electrolyte an ionic 
conductor by F" Frenkel type defects. The positive electrode is the 
electrode where the chemical potential of F 2 (g) is higher and the 



chemical potential of carbon is lower. This implies that 
3 

A/i c = — -A/^ Fr The cell reaction can be written as: 



Z - CrF 2 + Chopg + 3e = 3F- + i Cr 3 C 2 : ( + ) 

3F- + I Cr 3 C 2 = 3e + ? CrF 2 + C OL c : ( - ) 
Chopg = Cqlc 



(B) 



The main advantage of this kind of cells is that no supplementary 
thermodynamic data are necessary to derive the thermodynamics 
of the cell reaction. Only the emf vs T experimental data are 
necessary to obtain the reaction Gibbs free energy change, A r G. 
The reaction changes of enthalpy, A r H, and entropy, A r S, can be 
also obtained as shown in the followings. 

The adopted experimental procedure requires accumulating 
many isothermal and stationary values of emf through several stair 
shaped thermal cycles each one composed of tens of isotherms. 
Figure 1 shows the typical behaviour, almost equal for each one, of 
five tested cells. The first thermal cycles are characterized by 
hysteresis that tends to disappear. The presence of a vanishing 
hysteresis in the emf trend demonstrates that the whole system of 
the carbon shells in the OLC and OLC themselves change to 
reach new stable configurations that convert reciprocally and 
reversibly as function of temperature. These findings are 
qualitatively and quantitatively similar to the behaviour we found 
in the bundles of SWCNTs [1]. The main difference with the 
present experiment is the final shape of the emf vs T curve that is 
characterized here by an exponential growth-like behaviour. 

The trend of the Gibbs free energy change is given by: 



A r G T =-3J r E(T) 



(1) 



being 3 the number of moles of electrons exchanged in reaction B 
and T the Faraday constant. The curve E( T) is well fit by the 3 rd 
degree polynomial: 



E{T) = a + bT + cT 2 + dT 3 



(2) 



.The coefficients and related errors of fit of emf vs T are given in 
Table 1. The representation of eq 1 is given in Fig. 2 where at 
920.6 K the equality of the chemical potential of carbon 
^ c (HOPG) = ^ c (OLC)is satisfied. 

If each point of the curve of Fig. 2 were a point taken at 
constant temperature and pressure, the enthalpy change of the 
transformation, A r Hp, would be given by: 



A r H T - 



E(T)- 



(dE 
\df /F 



(3) 



The temperature dependency of eq 3 should be written as given in 
any chemical thermodynamics textbook by equation: 



8A r H T 
dT 



-ITT 



8 2 E 



--Ac P (T) 



(4) 



being Acp(T) = cp(OLC) — c>(HOPG) and Cp the molar heat 
capacity at constant pressure. In the case of reaction B and at the 
experimental temperatures, it is expected that Acp(T) should be 
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Figure 1. The electromotive force (em/) vs temperature of galvanic cell A. The progressive disappearing hysteresis is shown from 1 st to 3 rd 
thermal cycle (see §1.4 of the Supporting Information S1). The bottom plot combines the 3 rd and 4 th cycles where the red curve is the best fitting 
curve given by a 3 rd degree polynomial (see Table 1). The residuals of fit are also reported. Most of the error bars are within the size of the 
experimental points. 
doi:1 0.1 371 /journal.pone.01 05788.g001 



quite low. Equation 4 implies the same value of pressure at each should be considered where, 



temperature for both the electrodes, i.e., 



(8A r H T \ _dA r H T 
\ 8f )p~ dT ' 
The 1 st partial derivative against T of eq 3, i.e. eq 4, shows too 
high Acv(T) values which are unexpected for the process under 
study. For this reason, the assumption of constant pressure is not 
more valid. Therefore, the total derivative of E(T), i.e., 



dE _ /8E\ (8E\ dP 
dT~\If) p + \Jp) T df 



(5) 



cE 
IP 



1 (8A r G T 
3F\ 8P 



3J- 



-.AV 



(6) 



being AV= K(OLC) - F(HOPG) and V is the molar volume. 
By substituting eq 6 in eq 5 and combining with eq 3, the P, T 



dependency of A r Hj = 
found: 



-3T^E(T,P) 



dE AV dP 
df + 3Tdf 



]} 
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Table 1. Coefficients and related 


errors of emfvs T polynomial fit: E(T)=a + bT + cT 2 + dT i . 






alV 


6/VK 1 c/VK 2 


cffVK 3 


coefficient 1.207 


0.0045698 -5.7832x10~ 6 


2.437x10"' 


error 0.12433 


0.0004144 4.5905 x10~ 7 


1.691 x10~ 1 ° 


R = 0.99836; y_ 2 = 1 .5935 x 1 0 5 . 

doi:1 0.1 371 /journal.pone.01 05788.t001 



A,.H T 



3T{ E(T,P)-T 



dE AVdP 
df + 3Tdf 



(7) 



same quite high temperature where Cy is for both ones close to -R. 

Therefore, eq 12 becomes a quadratic equation in the 1 st 
derivative with respect to T of the pressure: 



Accordingly, eq 7 differs from eq 3 for the extra term 
AVdP\ dE (dE\ 

T Wdf)' which ent ^df^{df) P 
The 1 st derivative of eq 7 with respect to T is: 



dA,H 
dT 



--ITT 



d 2 E 
df 2 



dAVdP A dP , d 2 P 
T l!Tdf +AV df + TAV dT 2 - < 8 > 



From the total differential of AV(T,P), the relationship below 
follows: 



dAV 
dT 



8AV 



dAV 
~8P 



dP 

,df 



■■A(aV)- 



dP 
dT 



A(kV) (9) 



where ot and K are, respectively, the isobaric volume expansion 
coefficient and isothermal compressibility. By making use of eq 9, 
eq 8 can be conveniently modified in the form: 



dP 
dT 



t A(ocV)dP 1 
''A( K V)df~ A(kV) 



IT 



J 2 E 
df 2 



a, 2 V 



= 0 (13) 



The solution of eq 13 has been obtained by using the experimental 
d 2 E 

data, ~JY2' an< ^ ^ e knowledge of the T and P functions of the 

physical data of HOPG and OLC. The T functions are reported in 
Table 2 and their related coefficients given in Table 3. Unfortu- 
nately, a and K data of OLC are unavailable and we utilized, as 
first approximation, the data of fullerene Cgo [18]. Due to lack of 
data on the temperature dependency of the OLC molar volume, 
an estimation can be done starting from the definition of the 

isobaric expansion coefficient 0t,(T) = — y^^j ana ~ lls ^ 

dependency reported in Table 2 together the coefficients for C 6 o 
given in Table 3. Therefore, equation below has been considered: 



d& r H_ d 2 E f dP 
dT dT 2 \ dT 



AV- 



dP 



TAV- 



A(aV)- 
d 2 P 



dP 
df 



A(kV) 



(10) 



dT dT 2 

Equation 10 should be compared with well-known equation [17]: 



dA,.H 
dT 



--Acp(T) + [AV-TA(otV)} 



dT 



(11) 



which states that if some phase transition occurs, the equilibrium 
pressure does not remain constant as the temperature is varied. 

Since we can write that Acp(T) = Acy(T) + TA[ I, by 



equating eq 10 and eq 11, the following 2" order differential 
equation for P is found: 



AV 



dT 2 



dP 



dP 



Ac v 

~f~ 



-A(kK) - + 2A(*V) dT 



a 1 V 



-7,T 



dT 2 



(12) 



V(T)=V(T 0 )exp 



(do + d\ T)dT 



(14) 



L(r 0 )exp{ Ar 



rfo+y(r+r 0 ) 



V(To) 



where [T 



■ 7q) = AT(T+ Tg) holds. The final approximated 



value is justified by the amount in braces, which is = 3x1 0 -2 being 
AT= (T — To) <700K. The values of do and d\ are reported in 
Table 3. The V{Tq) value was taken equal to 7.31 xlO -1 m 
mol" 1 at T 0 = 298 K which is the V° value of OLC [19]) (see 
Table 3). Due to the lack of data on the pressure dependence of a, 
K and V° for OLC, all these parameters were considered constant 
with pressure for the carbon species present in both electrodes. 
This choice was mandatory to ensure a balanced behaviour to 

both OLC and HOPG. Figure 3 shows the trend of vs T 

in the left axis calculated according to eq 1 1 after solving eq 1 3 to 
dP 

find being the latter quantity plotted on the right axis. It 



dT 



dP 



should be noticed that the shape of the curve might be affected 
d 2 P 

by the approximation — j =0 introduced in equation 12. 



To simplify eq 12, we can reasonably neglect an d assume 
Acy = 0 being both OLC and HOPG graphitic materials at the 
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Figure 3. Representation of — vs 7~obtained by solving eq 13 (right axis) and trend of A,.H(P,T) calculated by eq 11 (left axis). 
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From the total differential of A r H(T,P) the quantity 
dA,H fBArlA 
fdA r H\ dT ~ V dT 



\ 8P 



dP 
dT 



can be obtained and plotted 



fdA,.H T \ 

against T in Fig. 4 (right axis) together I — — — J (left axis). 

Equation 4 was used for this calculation. 

The entropy change of transformation B as function of 
temperature and pressure, was calculated as: 



A r S{ T,P) = [A r H( T,P) —A r G( T,P)]/ T 



(15) 



.This equation can be calculated through eq 7, after solving the 
quadratic eq 13, by using the experimental data of A r G reported 



on Fig. 2. These thermodynamic functions are shown in Fig. 5. At 
920.6 K, A r H and A r S values are 95.8 kj mol" 1 and 104.1 JK" 1 
mol -1 , respectively. 

Why does the pressure change with temperature in our 
experiment? The experiment is carried out under high vacuum 
and the static pressure on the cell is maintained constant at about 
0.1 bar as described in the experimental section (see § 1.3 of the 
Supporting Information SI). Thus, the present experimental 

dP 

conditions cannot justify the presence of the — quantity. The 
dP 

explanation for > > 0 should be searched in the growth of the 

pressure inside the cell, which is caused by a process occurring 
under volume constraint. The positive sign of emf indicates that 
|((OLC) > /((HOPG) where [X is the chemical potential of carbon. 
By itself this does not justify the increase of P with T but the trend 
of emf with T suggests it. Furthermore, there are experimental 



Table 2. Temperature functions of the molar volume, V°, isobaric volume expansion coefficient, a, and isothermal compressibility, 
k, used for computing eq 13. 









«(7) Ibl 




HOPG 25 


a\ +CI2T + CI3T 2 


a 2 +2a^T 


(b l +b 2 T+b 3 T 2 y' 






a\ +CI2T + CI3T 2 


OLC 


V° 







[a] Fullerene C 60 data were assumed for OLC from reference [18] (see also text) with the exception of which is known [19]. 
[bl Linear equations for OLC (C 60 ) were taken up from the plots in reference [18]. 
doi:1 0.1 371 /journal.pone.01 05788.t002 
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Table 3. Values of the coefficients of Table 2. 



Coefficient Ial 


Value 


Ref. 


i?i/m 3 nriol 1 


D.oU X 1 u 


LzbJ 


^2' m mol l\ 


z. 1 4 X 1 U 


[25] 


(3 3 /rn mol i\ 


1.95x10 14 


[25] 


bi/Pa 


3.63x10 10 


[25] 


bJPa K 1 


-1.32X10 7 


[25] 


&,/Pa K 2 


9.45 xlO 3 


[25] 


to/Pa" 1 


8.37x10"" 


[18] 


c,/Pa _1 K 1 


7.66x10~' 5 


[18] 




4.30 x10~ 5 


[18] 


£/,/K 2 


4.27 x10~ 9 


[18] 


l/Wmor 1 


7.31 x10~ 6 


[19] 



[a] The original values of coefficients a, are given in nm 3 . They were multiplied by 10 27 W A being N A the Avogadro number. 
doi:10.1371/joumal.pone.0105788.t003 
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evidences indicating that the electrode containing OLC is the site 
where the volume constraint is particularly high. They arise from 
the comparative analysis of both electrode powder mixtures, 
before and after the experiments, by making use of x-ray 
diffraction (XRD), micro Raman spectroscopy (mR) and high- 
resolution transmission electron microscopy (HR-TEM): 

• The XRD patterns spectra demonstrate that no chemical 
change occurred (see §1.5.1 of the Supporting Information 
SI). The large increase of crystallinity of CrF 2 was the unique 
difference found in XRD. Due to their almost amorphous 
state, OLC give low intensity XRD features differently from 
HOPG, which is very well revealed. In the present case, the 
situation is complicated by the presence of features of other 
species in the electrode that superimpose and/ or they are very 
close to the OLC features (see §1.1.2.2 of the Supporting 
Information SI). It is reported [20] in high temperature - high 
pressure experiments on OLC, that their most intense feature 
0200 = 0.354 nm, 20 = 24.73°; see §1.1.2.2 of the Supporting 
Information SI) is negligibly shifted as function of pressure at 
constant T. For instance, at 500 °C the lattice parameter 
changes with pressure as —9.6x10" ' nmPa . At our 
calculated pressures, which are in the order of tens of MPa 
(see below), a shift of = 10~ 5 nm is expected, which is clearly 
undetectable. The mR spectroscopy does not reveal differences 
(see § 1.5.2 of the Supporting Information SI) with exception of 
some changes in the relative intensities, which are meaningless. 
There is superimposition of the OLC and HOPG features with 
the features of the mixture Cr 3 C 2 - CrF 2 (see §1.1.2.3 of the 
Supporting Information SI). Therefore, though this inconve- 



nient, the absence of any new feature is a proof that both 
carbon species remain unchanged; 
• HR-TEM and STEM analyses show a very different scenario 
of the electrode powder with OLC before and after the 
experiments. Before, OLC are easy identifiable in the mixture 
as reported in the STEM image of Fig. 6A. Free OLC 
nanostructures (red arrows) are quite abundant, where they 
can be seen either as spots with very light contrast on the holey 
carbon grid or mixed with the Cr F 2 powder. The 

HAADF-STEM mode is also known as Z contrast mode 
because the intensity in each point of the image is proportional 
to Z" with a~2 and Z the atomic number. Carbon appears 
dark, as it is much lighter than the other compounds present in 
these samples. The OLC morphology at the nanoscale is 
shown in Fig. 6B and is found in agreement with previous 
experimental results [7] [21]. The inset in Fig. 6B shows the 
Fast Fourier Transform (FFT) of OLC where the reflections of 
graphitic planes appear quite diffuse. 

The powder mixture of the same electrode after the experiment 
appears completely different: OLC are practically disappeared 
while large grains with round shaped contours were found as 
shown in the STEM image in Fig. 7A. Free OLC are hardly 
observable: just small pieces on the agglomerates, while the holey 
carbon grid is practically clean. The analysis on the very few OLC 
found in the sample reveals some significant differences with 
respect to the OLC before the experiment: the graphitic planes 
form now better ordered quasi-spherical structures (see Fig.7B). 
The electron energy loss spectroscopy measurements (see § 1.5.3. 
of the Supporting Information SI) agree with this result: the sp z 
carbon coordination values were found 82 ±3% and 86±3% for 
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Figure 6. Before the experiment. Panel A: STEM image of the powder of the electrode containing OLC. The OLC nanostructures appear on the 
holey carbon grid as spot with very light contrast. Some of them are indicated with red arrows. Panel B: HR-TEM of OLC. The FFT image in the inset 
shows poor crystallinity of OLC. 
doi:1 0.1 371 /journal.pone.01 05788.g006 
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Figure 7. After the experiment. Panel A: STEM image of the powder of the electrode containing OLC. The OLC nanostructures are practically 
disappeared. The holey carbon grid is clean and only large grains with rounded shape stay on the grid. A few visible OLC are indicated with red arrow. 
Panel B: HR-TEM of OLC. The FFT image in the inset shows that the crystallinity of OLC is higher with respect to the FFT image in Figure 6B. 
doi:1 0.1 371 /journal.pone.01 05788.g007 
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Figure 8. Zero loss TEM image (Panel a) of a typical small sintered particle (as the particles in Figure 7A) of the electrode with OLC 
after the experiment. The particle contains several cavities, that can be empty, or at least full of material much lighter than the surrounding 
sample, as it is proven by the comparison between the HRTEM (panel b) and STEM (panel c) images of the same cavity, that is whiter in b and darker 
in c. 

doi:1 0.1 371 /journal.pone.01 05788.g008 



OLC before and after, respectively and the result seems to point to 
a further graphitization of OLC towards multishell ordered 
structures in the process (see discussion below). The sintered 
agglomerates usually are so thick to result non-transparent to 
electrons and to prevent the analysis of their internal structure. 
When small particles are found, the EELS revealed typical features 
of variable size inside the particle. Their contrast is typical of 
cavities, but it was not possible to assess whether they are empty or 
possibly filled with OLC. This is the most reasonable hypothesis, 
being OLC otherwise missing. 

The zero loss TEM image of a typical small sintered particle 
containing several cavities is shown in Fig. 8a. The fact that these 
cavities are empty or filled with material lighter than the 



surrounding sample is proven by the comparison between the 
HRTEM and STEM images in Fig. 8b and Fig. 8c of the same 
cavity, that is whiter in b and darker in c. 

The electrode containing HOPG and the same Cr 3 C2 - CrF 2 
mixture does not show any modification after the experiment, as it 
can be seen by comparing the STEM images in Fig. 9 A and B. 

The following considerations arise from the above evidences: 

• The temperature range explored in the experiments does not 
allow suspecting any vaporization of carbon from OLC. The 
C-nanostructures are still in the electrode; 

• High temperatures and experiments lasting hundreds of hours 
create the ideal conditions to favour the sintering of the whole 
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Figure 9. STEM images of the electrode containing HOPG. Panel A, before the experiment. Panel B, after the experiment. No differences 
were found and no sintering was observed. 
doi:1 0.1 371 /journal.pone.01 05788.g009 



mixture where OLC remain embedded in the very hard and 
refractory structure of Cr 3 C 2 . The melting points of Cr 3 C 2 and 
CrF 2 are respectively 2168 and 1167 K. The highest 
temperature reached in the experiments was 1030 K, which 
is close to the melting point of CrF 2 Thus, it is expected that its 
softening acts as a sintering medium of the whole mixture. The 
high specific surface area of OLC (>100mg , see § 1.1.1 of 
the Supporting Information SI) contributes to this process and 
facilitates the embedding of the nanostructure. Therefore, the 
role of OLC nanoparticles seems to be crucial and it explains 
the absence of any modification of the electrode containing 
HOPG (and the same Cr 3 C 2 -CrF 2 mixture) which worked in 
the same experimental conditions; 

Consider as rough approximation the electrodes as cylinders 
free of expanding with temperature filled with rigid particles 
(sintered Cr 3 C 2 -CrF 2 mixture) embedding nanoparticles 
(OLC) the dilatation of which is hindered. By imposing no 
volume change for both OLC and HOPG, the respective 
dP a 

— — = — values can be evaluated and compared as follows: 
dT k 



dP 
dT 



-) s0.5MPaK 

kJ o 



.[ dP 
'df 



-) S^PaKT 1 (16) 

H VK/H 



The subscripts O and H stand for OLC and HOPG. The 
temperature dependence of the quantities entering eq 1 6 has been 
neglected in the computation and the coefficients ffl 2 , b\, c 0 and do 
were only retained in the expression reported in Table 3. The 
values found in eq 16 demonstrate clearly the great difference in 
the internal pressures to which both carbon forms are subjected in 
the respective electrodes. For the sake of comparison, the quantity 



dP 



dT 



for Cr 3 C 2 is: 



dP 
df 



s3.4MPaK _1 where k [221 and 

Cr VK/ Cr 

a [23] are available in literature. Due to the negligible changes 
with temperature of the internal pressure (see eq 16), the chemical 
potential of carbon in HOPG, /i H , is only a function of T 
differently to the chemical potential of carbon in OLC, n Q , that is 
a function of T and P(T). Therefore, we can write that: 



A r G(T,P)=no[T,P(T)]-n n (T). 



(17) 



Evaluation of the surface contribution 

Owing to the nanometric size of the carbon particles, the 
surface contribution to the energetics of the process cannot be 
neglected, in principle. In the following of this section the shape of 
the nanoparticle is assumed to be spherical. Let us estimate the 
chemical potential difference using the thermodynamic approach, 
based on the Gibbs model for interface [24], in terms of bulk 
chemical potential, surface free energy and surface tension. 

The difference between surface energy and surface tension plays 
a fundamental role when dealing with interphases in the solid 
state. In this respect, we briefly recall the definition of these 
quantities, which are needed for the discussion that follows. The 
surface energy excess, a, is the reversible work required to create a 
unit surface at constant volume, temperature, pressure and 
composition [24]. In such a process the surface is not stretched 
during its formation. On the other hand, the surface tension, y, is 
defined as the reversible work (per unit area) required to change 
the area of the surface through a stretching process. It is worth 
noticing that in this definition the infinitesimal deformation of the 
surface involves an "initial" configuration of the system where the 
solid is not stressed. By including bulk and surface contributions, 
the general expression of the chemical potential difference reads 
(see §2 of the Supporting Information SI) 



A,-G( T) = A/i B + Afi s = A/i B 



2 V 



ky 
3r 



(18) 



where r is the radius of the nanoparticle and A/( B denotes the 
"bulk" contribution, that is A^,.^^. In the following we focus our 
attention on the impact of the surface term (A/i s ) on A r G(T) and 
on its temperature dependence. 

The contribution of the surface tension containing term to Afi s 
is usually smaller than that due to the surface excess free energy, 
2Va 

. Consequently, A/i s is expected to be higher than zero. 

Typical values of the quantities entering eq 18 are [25] 
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K=5.36xl0~ 6 m mol" 1 , K = 3.24xl0" 10 Pa" 1 and r= 10 nm. 
As far as the a value is concerned, for nanocarbon structures, such 
as fullerene, single and multi walled nanotubes, it is found [5] to 
range between 0 and 0.045 Jm 2 . On this basis, the upper bound 

2K<7. . _! 

of Ap s = is estimated to be about 50 Jmol for r= 10 nm. 

As far as the temperature dependence of Afi s is concerned, it is 
ascribed to the temperature dependence of the physical quantities 
entering in eq 18. In particular, the temperature dependence of <x 
is given by the Gibbs adsorption equation that, for one component 

system, implies = — s* , where s x = S x / A is the excess surface 

entropy. From this equation, neglecting as a first approximation 
the temperature dependence of the excess entropy in the 
temperature range here considered, we get 
a(T) = —s x (T—T*) + (j* where T* is the temperature at which 
A r G(T)=0. Concerning the temperature dependence of the 
surface tension, given by eq 1 in the Supporting Information S 1 (§ 
2), it is obtained 



dy 
df 



-s x - 



J d s s x 2 
d.A 3 



(19) 



where we exploit the fact that the surface energy is a function of 
the independent variables temperature and "stretched area" A (the 
subscript s stands just to remind us that the variation is performed 
by stretching). According to the discussion that follows this 
derivative can be estimated using the Born-Stern method [26] . In 
this approach, the entropy excess per unitary area of a solid- 



vacuum interface is equal to s" = — (^S — S j = p a [s\ ~ s \*j where 

S is the entropy of the actual solid, S' the entropy of the model 

system (without interface), p a the surface density and S\ (s^j the 

entropy, per atom, of surface (bulk) species. By assuming, as usual, 

^Si ~ s i^ independent of temperature and deformation of the 

F d s A d s p a p a 

surface — — , lor 



d s A A 
d s p a 2 2 



one obtains 



x J ' \ "*Pa 1 1 

■■ —s — A \s\ — Si 1 -j— + -txy= - ay. For instance, by i 



dy 
df 

sidering the excess entropy equal to the change in the vibrational 

v± 

entropy of the atom the expression s* = kp a In — holds where 

\ 

v^(vj_) is the vibrational frequency of the surface (fully 
coordinated) atom along the "broken bond" direction. Since for 
the number of mole of atoms in the particle, N, is constant we set 



— = -i-N 
V 4n 

Moreover, 



= b 1 in eq 18 to obtain Ap s = 2br 2 - 



Kf 



dAp s 

~dT 



r 



a /„ 5 kv 2 



(20) 



where a is of the order of 10 6 K '.A crude estimate of the 
vibrational contribution of S x , on the basis of the Kossell model 

r 2 

[26] shows that s x >kp a In v2 that implies, using p a = pi with p g 
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Figure 10. The change of internal energy AU[T,P(T)] of OLC is reported against temperature according to eq 22. 

doi:10.1371/journal.pone.0105788.g010 
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the graphite density, a value of s x > 10 4 Jm 2 K 1 . This figure is 
about two orders of magnitude larger than oc. Moreover, 
(2K?7r) r=10nm >0.1 J mol" 1 K" 1 . From eq 15, the value of 
A r S(P,T), in the middle of temperature range, is equal to 
96.75 Jmol _1 K _1 , which is not consistent with the value 
(-2Vs x /r). 

The present analysis indicates that the bulk term is the 
dominant contribution both to the chemical potential change of 
transformation B and its derivative on T. At this point, a comment 
is in order on the use of the Gibbs model for dealing with OLC 
systems. In fact, in this approach "bulk" and surface terms sum up 
as two independent contributions where the "size effect" is usually 
contained in the surface term only, i.e., Aji B = A[i r ^ rxi is 
independent of r. On the other hand - still remaining in the 
framework of the Gibbs model - owing to the variable curvature of 
the sheets, which made up the nanoparticles, in the OLC the 
"bulk" term has to be considered a function of the particle radius. 
For this reason in the following we refer to "atomistic" approaches 
for determining the energetics of the OLC particles. In fact, in 
these methods a suitable interaction potential functions for the C 
atoms is used and this makes it possible to estimate A/i B as a 
function of r. 

Now we focus our attention on the internal energy of the 
particle due to the isothermal transformation from a relaxed state 
at internal energy U u and volume Vq to a stressed state at U and V. 
To do this the energy equation, namely 



8U 
dV 



rr- 



was utilized. The integration of eq 2 1 gives the change of internal 
energy as (see §3 of the Supporting Information SI): 



AU= (U- U Q ) T = — { [1 - (1 + KP)e- KP ] + ocT(e- KF - l) } (22) 

K, 

The value of v, a. and jck are given in Tables 2 and 3. The 
pressure as function of T was obtained by numerical integration of 
dP 

= f(T) assuming P = 0 at T= 776.6 K, which was the lowest 

experimental temperature. The resulting function satisfies the 
relationship P = mo+m\T + m2T 2 +m^T 3 where 

m 0 = 1.2947xl0 6 bar; w*! = —4187.2 bar K" 1 ; m 2 = 4.2571 bar 
K -2 ; ?»3= — 1.3067xl0 -3 bar K~ 3 with correlation coefficient 
0.9999. At the highest temperature, the maximum pressure value 
is about 70 kbar. The trend of AU[T ,P(T)\ with temperature, 
which is almost perfectly quadratic, is plotted in Fig. 10. 

It is worth noticing that the AU[T ,P(T)\ values and its trend 
with temperature are not comparable with values and trend of 
A r H(P,T) reported in Fig. 5. For example, at 903 K, which is the 
temperature in the middle of the explored temperature interval, 
A r H(P,T) 



the ratio 



AU[T,P{T)] 



■ 10 



0.30 



0.20 



0.10 



0.00 




750 



Figure 1 1 . Trend of the change of the fraction of C defects in OLC against T. One assumes that the real transformation involves equilibrium 

C H opg = Colc, =J'C fl rf + CoLc v _, and y is calculated by eq 23. 
doi:1 0.1 371 /journal.pone.01 05788.g01 1 
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Due to such a high ratio jointly to the very small surface 
contribution, we assume that as a consequence of the high internal 
pressure another process connected to B should be considered. We 
are likely dealing with a reversible process that requires the 
formation of C defects in OLC, probably already defective. Thus, 
the whole process determining the emf of the cell could be written 
as: 

Chopg = C 0 lc x =yCad J rCo\ J c x - y (Bl) 

instead of process B. 

Obviously, this implies that the previously considered relation 
A r G = fi 0 — should be written as A r G= < fi 0 > — n H where 
<)1q> is the mean chemical potential of carbon atoms in the 
electrode containing the OLC nanostructures. In this quantity, the 
C atoms formed should be considered adsorbed atoms, C a d the 
formation of which causes dangling bonds. A rough evaluation of 
the defect fraction in equilibrium Bl can be evaluated through 
equation: 

= \H(P,T)-\ f H 0x _ y ^ A r H(P,T) m) 
y " AfH c ' ~ A/He [ ' 

where AfHc and AfHc are the enthalpy of formation per carbon 
atom of defective OLC and free carbon atom (7.43 eV at 298 K) 
[27], respectively. The quantity AfHo x _ y depends on the number 
of shells as described for multishell fullerene (MSF) made of 
concentric circumspheres of polyhedra C„ (//,, n = 60k 2 ) and C n (//,, 
n=180A 2 ) where k =1,2 up to L is any positive integer 
representing the number of shells, lu and n are respectively the 
symmetry group and number of C atoms. According to the 
literature [28], the quantity AfHo x _ y is of the order of magnitude 
of meV and therefore it is negligible when compared to A r H(P,T) 
which is about 1 eV. Even the formation energy of a single carbon 
shell starting from graphene implies negligible energy values 
ranging from 331 to 34 meV at 1 for C 60 and C 960 , respectively 
(see J. Bernholc et al. [3]). As stated in the above reference, the 
authors showed that the larger the size of the shell, the lower the 
energy of formation. Figure 1 1 shows how y changes with T. 

Due to the assumption for obtaining eq 23, the trend of y is the 
same as A r H(P,T) shown in Fig. 5. Anyway, the value of y is 
likely consistent with a change corresponding to one or more OLC 
shells. The experimental evidence through the emf vs T data shows 
that this process occurs reversibly. In addition, defects imply the 



formation of dangling bonds, which also contribute to the increase 
of the entropy change. 

Conclusions 

In the present experimental work, we obtained high tempera- 
ture data of the stability of onion-like carbon with respect to highly 
oriented pyrolytic graphite. Since this was performed in reversible 
way, through the emf measurement of a galvanic cell, the 
quantities related to the transformation investigated have to be 
considered reliable thermodynamic data. 

The evidence that in the operating conditions the emf data 
cannot be only function of temperature allowed measuring the 
changes of the chemical potential of carbon in OLC under volume 
constraint. The volume constraint is reasonably proven by the 
HR-TEM images where the rigid cages made of sintered CrF 2 - 
Cr :i C 2 mixture likely embed OLC. 

The positive value of A r S and its increase with temperature 
indicate that the OLC is a nanostructured system with a large 
number of defects. The change of their fraction with T was 
evaluated. The present analysis indicates that the bulk term is the 
dominant contribution to both chemical potential and entropy 
change of the transformation. 

The absence of any chemical change in the electrodes after 
hundreds hours of work at high temperature guarantees the whole 
reliability of the experiments. 
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